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Heavy baryon chiral perturbation theory for pentaquarks is applied beyond leading order. The 
mass splitting in the pentaquark anti-decuplet is calculated up to NNLO. An expansion in the cou- 
L^ ', pling of pentaquarks to non-exotic baryons simplifies calculations and makes the pentaquark masses 

r— ^ ' insensitive to the pentaquark-nucleon mass difference. The possibility of determining coupling con- 

("—^ ' stants in the chiral Lagrangian on the lattice is discussed. Both positive and negative parities are 

fvj , considered. 

^ ' I. INTRODUCTION 

lO ■ Recent experiments have found evidence for an exotic hadronic resonance, the 0+(154O) [u IM 13. 13. la. la ItL Isl 19. [lOl| . 

This has initiated an effort in identifying other members of the same multiplet and answering a host of questions 
^^J \ about these resonances. What are their masses and widths? Why are their widths so small? What is their spin and 
J^ ' parity? How are they produced? How certain is their existence? 

«y~, ' Various models have been proposed for the pentaquarks that answer these questions with varying degrees of success. 

fvq It would also be helpful to have a set of model independent results using a minimal set of assumptions common to 

\^ I the models. This suggests developing an effective field theory for pentaquarks. Since processes of interest, such as 
^^ 1 pentaquark decays, produce nucleons and soft mesons, it is natural to consider chiral perturbation theory including 
■^ ' these fields. 

^^ In '22|, chiral perturbation theory is applied at leading order to decays of pentaquarks. However, in the formulation 

r^ ' of the theory used, baryon masses in the Lagrangian make the power counting beyond leading order difficult. With 
O^i some tweaking, namely the treatment of baryons as heavy fermions, the chiral perturbation theory offers a simple 
JL. ■ power counting and ease of calculation of Feynman diagrams |l9l | . In this form, it has been applied at leading order 
fl^ ' to decays of pentaquarks of both positive and negative parity and of both spin ^ and | |23l | . The results for positive 
,-^ ■ parity and spin i agree with those in |23| up to higher order corrections, as expected. 
^ I In the past, heavy baryon chiral perturbation theory has been applied to the octet and decuplet of baryons and 

•rH . to heavy mesons. The theory has been used to systematically calculate next-to-leading order (NLO) corrections in 
^% ■ nucleon-nucleon scattering that agree with experimen t IlIL ri2| . The theory has also been used to calculate corrections 
to the Gell-Mann-Okubo formula for baryon masses |23|.The results have been useful in the chiral extrapolation of 
lattice predictions for baryon masses [13]. 

In this paper, we apply heavy baryon chiral perturbation theory to pentaquarks beyond leading order. We calculate 
the mass splitting in the pentaquark anti-decuplet up to NNLO. We also identify an expansion in the coupling of 
pentaquarks to non-exotic baryons that simplifies calculations. 

Throughout, we assume that 0"'"(154O) is a member of an anti-decuplet. Experiments suggest that 8+(1540) is an 
iso-singlet with hypercharge two. The smallest flavor multiplet that can accommodate it is an anti-decuplet. Such a 
multiplet of pentaquarks appears naturally in the diquark model, along with a degenerate octet of pentaquarks [T^ . 
An anti-decuplet also appears in chiral soliton models ^g and in uncorrelated quark models, along with several other 
multiplcts of pentaquarks of different masses. In this paper, we ignore the possibility of other pentaquark multiplets 
with a similar mass. We will take up the consequences of including other pentaquark multiplets in a later work. For 
most of the paper we assume that the pentaquark fields have spin one-half and positive parity. However, we also state 
the results for negative parity to show that the mass splitting does not distinguish between parities. 

In section ^ we review heavy baryon chiral perturbation theory for pentaquarks and develop a small coupling 
expansion. In section llTTl we develop the power counting for this theory and determine the diagrams contributing to 
the self-energy at NNLO. In section llVl we calculate the mass splitting in the pentaquark anti-decuplet at NNLO. In 
section IVI we treat negative parity. The final section has conclusions and directions for future work. 
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II. CHIRAL LAGRANGIAN AND SMALL COUPLING EXPANSION 



The usual chiral perturbation theory for the pseudo-Goldstone bosons can be supplemented by the CCZW method 
to incorporate other "matter fields" into the theory |17l [la | . The method prescribes the most general Lagrangian for 
matter fields, that is invariant under chiral symmetry. We briefly review our field definitions and conventions and 
then construct the desired Lagrangian. The pseudo-Goldstone bosons are described by the field ^ = exp(i7r//) where 
/ = 131 MeV is the pion decay constant and the field it is 
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Under a chiral symmetry transformation, the field ^ transforms as 

^{x) -^ L ^{x) U\x) = U{x) C{x) R\ 

where U{x) defined in terms of L, R, and £,{x) by eq. ^. 

The octet of baryon fields, which we refer to as nucleons from now on, is 



B = 



Under a chiral symmetry transformation, the field B transforms as 

B{x) -^ U{x)B{x)U{x)'' 
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where U{x) is the transformation that appeared in eq. 1211. In fact, the CCZW results show that any transformation 
for B that reduces to the adjoint transformation for the unbroken vector subgroup is acceptable. The different choices 
simply correspond to redefinitions of the baryon field. 

An anti-decuplet representation SU(3) is the completely symmetric part of the tensor product of three anti- 
fundamentals. The components of the pentaquark anti-decuplet, with respect to the quark basis, are organized 
below according to the weight diagram. 
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We will drop the 10 subscripts when there is no possibility for confusion. Under chiral symmetry, the field P transforms 



as 



P^Jk{x) ^ Pl.mn{x)U\xy,U\xr,U\xrk 



(6) 



To construct a chiral Lagrangian for matter fields, it is convenient to trade in the ^ field and its derivative for two 
vector fields: 



A, = \{^d,e~ed,0^jd,^~^^ 



3 ild^Ti,-!:] ,7r] 
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These vector fields have the simple transformation properties: 

A^{x)^U(x)A^{x)U^{x) 

V^{x) -^ Uix)V^{x)U\x) - d^UU\x). 
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In QCD, the quark mass Mq explicitly breaks chiral symmetry. However, if Mq is assumed to transform as 

Mq -> LMqK^ (11) 

then the QCD Lagrangian is invariant under chiral symmetry. In order to have the same symmetry breaking pattern 
in chiral perturbation theory, we have Mq transform in the same way. It is also convenient to trade in Mq for a 
combination of Mq and ^ which transforms like the other fields: 

Ml = {£}Mq£}+iMl^) (12) 

M| ^ U{x)mIu\x). (13) 

A covariant derivative can be defined for nuclcon and pentaquark fields by 

D,,B^d^B+[V^,B] (14) 

{D^P)^^^ = d^P,,k - P^Vl, - PukVl, - P^.iVlu (15) 

As usual, the covariant derivatives of the fields transform like the fields themselves. Under parity, the fields B, V^, 
and M^ are even, while the field A^ is odd. Both even and odd parity will be treated for the field P. For most of the 
paper, we will assume that P is even under parity. The results for odd parity will be summarized in section IVl 

Having discussed the symmetries of the theory, we now focus on the relevant scales and the power counting in the 
theory. Contributions from terms in the chiral Lagrangian with derivatives acting on baryon fields are suppressed by 
?71b/A^; however, this is not a small quantity. It is necessary to sum these contributions to all orders in tti^/A^. 
The difficulty of doing this in practice suggests using an alternative approach called heavy baryon chiral perturbation 
theory (HBxPT) jl^. It uses well-tested ideas from heavy quark effective theory. In HBxPT, the momenta of baryons 
are split up as 

p)" = msu'' + fc'', (16) 

where v ■ k <^ rriB in low-energy processes. In other words, the baryons essentially have constant velocity when 
interacting with soft mesons. The usual baryon field B{x) is traded in for baryon fields of definite velocity 

B„(x) = i(l + ^)e^™-'^-i3(x), (17) 

which satisfy ^B^ = By. This redefinition cancels the tub term in the leading order chiral Lagrangian resulting in a 
simple baryon propagator i/{v-k). Moreover, since the baryons are not far off-shell in low-energy processes, derivatives 
acting on By produce powers of the small residual momentum k, which is of the size of the meson momenta, rather 
than the total baryon momentum p. Thus, a simple power counting is restored, in terms of the small quantities k/A^ 
and k/niB- 

In pentaquark decays, there are two multiplets of matter fields with different masses. Applying this procedure for 
each field separately results in simple kinetic terms for both fields but introduces a phase in the interaction terms. 
An equivalent possibility is to use the same field redefinition for both fields, 

P„(x) = i(l + ^)e™^^-^P(2:) and B„(x) = ^ (1 + ?^) e""^''-^B(a;), (18) 

where mp is the mass of the pentaquarks, to keep the interaction terms simple. Doing so introduces a constant 
A = mp — TUB in the nuclcon propagator, i/{v ■ k + A). We adopt this second approach, but the results are the same 
using either approach. 

In terms of the fields Py and By, the lowest order chiral Lagrangian is £ = Cp + Cpg + Cg + C-„ ■ The usual 
pieces involving only nucleons and mesons are, 

6g^ = Tr {By {ivD + A) By) + 2D Tr (B.S^f {A^, By}) + 2F Tr (B„5^ [A^, By]) (19) 

£(2) ^ L Tr (^M^^) + Bo^ Tr Af |. (20) 

Here S^ is the (Pauli-Lubanski) spin operator in a frame with four-velocity v. The constants D and F are the usual 
axial vector coupling constants with gA = D + F. The pieces of the Lagrangian involving pentaquarks are 

C^p'> = PyiV ■ DPy + 2HpPySy ' APy (2l) 

L^pI = 2CPB (PvSy ■ ABy + BySy ■ APy) . (22) 



FIG. 1: a) A pentaquark-pion-pentaquark vertex, b) A pentaquark-pion-nucleon vertex. This is the only diagram contributing 
to pentaquark decay at leading order in the chiral expansion, c) A quark mass insertion. This is the only diagram contributing 
to the pentaquark self-energy at NLO. Double lines are pentaquarks, single lines are nucleons, and dashed lines are pions. 
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The constants CpB and Tip are pentaquark analogues of constants introduced in T!J| for the baryon decuplet. The 
contraction of flavor indices has been suppressed because there is only one possible contraction between the multiplets 
that results in a singlet. More explicitly, 



Cf - P'^Hv^' {D^P)^^^ + 2'HpP'''S^AliP,k (23) 

£% = 2CpB {P''''S>iA^\B^,ekirn + B^'A^ I'S^P.^kc"'"') (24) 



To be completely explicit, we should also denote the Lagrangian £ above by Cy. The full Lagrangian can then written 
as £fuii = / ■2^^v- Since we will be working in a single velocity sector, we drop this formality. Finally, note that C-n 

would take its more familiar form in terms of the field S = ^■^. 

The only constants with dimensions of mass appearing in this Lagrangian are Mq and A. In chiral perturbation 
theory, Mq is treated as second order in the chiral expansion because m^ ex ruq. The mass splitting in each multiplet 
of matter fields due to the quark masses happens at next to leading order. On the other hand, the mass difference A 
that appears in the Lagrangian does not vanish in the chiral limit; formally A oc A^. At NLO in chiral perturbation 
theory, the mass difference can be expressed as 

/ ( 2m\ + ml W ( ( 2m\ + ml W 

using eq. (|32|) below for the NLO Sjq mass and the well-known analogue for the S nucleon. To obtain a rough estimate, 
we use the mass of the candidate S seen in 6] at 1860 MeV and assume that <Tp and up are approximately equal. 
We find A -- 750 MeV. Throughout this paper, we take A^ = 47r/ = 1.65 GeV so that A/A^ -- 1/2. This is a 
borderline situation that could lead to problems. Fortunately, A does not appear in the mass splitting up to NNLO. 
Thus, there are no problems in this calculation. 

The partial widths for the decay of a pentaquark can be easily calculated from the Lagrangian at leading order. 
Only the tree diagram in figure^ contributes at this order. The amplitude is 

A{P'' ^ B^n') = F-^-^^B^fc^ • 5„P„ (26) 

where F"'^^ is a flavor factor. The partial width is then 

r(P'^ ^ B^^-^) - ^ / dlifSih ~kf){\Y.\AA (27) 

= {F^n'^,—kl (28) 

In evaluating the total width of 0+ numerically, we use the mass of the O for mp and the mass of the proton for iriB , 
taking into account some NLO corrections. The flavor factor [F"-'^'^)'^ for both decay channels of the 8+ is 1 and the 
momentum of the kaon produced is 269 MeV. The total width of 0+ is 

r(e+) = r(e+ -^ k^p) + r(e+ -^ K+n) = (220 MeV) clg. (29) 



There is one minor difference between this value and that in eq. (10) in [23 due to a choice for the value of / used 
in eq. 1)2 7|l . In ordinary chiral perturbation theory, one often finds that in processes involving kaons, some of the 
one-loop corrections can be accounted for by replacing / by Jk in the tree level result. Without calculating to higher 
order, it is not clear that this improves the accuracy for pentaquark decays. Thus, we use the leading order / for all 
decays. 

Experiments have found the total width of 0+ to be surprisingly small. In several experiments, the width is 
below the detector resolution. Two recent experiments suggest that the width is in the neighborhood of 10 MeV and 
possibly even smaller 8, 10]. Reanalysis of old high-statistics KN scattering data and the non-observation of 0+ 
there suggests an even smaller width, possibly smaller than 1 MeV |2J, [23 . The recent experimental results on the Q 
width suggest that CpB is at most j, and the reanalysis of old data suggests that CpB is smaller than -V. The small 
coupling constant Cps allows a further expansion. Thus, we now have an expansion in the small quantities fc^r/A^, 
k-nimp, and CpB- 

The small coupling argument, made for the nucleon-meson-pentaquark vertex, is also relevant for the coupling 
of pentaquarks to other observed baryons with masses near that of pentaquarks. Based on dynamical assumptions 
coming from quark models, we expect non-exotic states with a baryon and meson to have small wavefunction overlap 
with a pentaquark |l5l l26l l27l|. For example, in the diquark model, one of the diquarks in a pentaquark has to be 
split up in order to create a non-exotic state. 

On the other hand, the pentaquark-meson-pentaquark coupling Tip should not be suppressed by this effect. This 
is in analogy with the couplings g^ and TL for the nuclcons and the dccuplct respectively, which have been fit to data 
and found to be of order 1. 

Models are needed at this stage to provide an estimate of coupling constants not measured in experiment. Chiral 
perturbation theory is definite about the fiavor, spin, color, spatial properties of hadrons. For example, flavor Clebsch- 
Gordon coefficients between hadron multiplets can be calculated directly. However, chiral perturbation theory says 
nothing about the wavefunctions of quarks in the hadron. That is the realm of QCD. In this way, QCD determines the 
coupling constants appearing in the chiral Lagrangian. Although QCD itself is remarkably simple, it is complicated to 
compute things, such as the coupling constants. Quark models attempt to approximate this complicated calculation. 

The small coupling argument is relevant for pentaquarks of negative parity also. The Lagrangian Cp is unchanged 
under a change in parity of the pentaquark field because the field appears quadratically in each term. The Lagrangian 
CpB, however, is affected. The chiral Lagrangian (before the heavy baryon transformations) changes by the replace- 
ment of 7'^7^ with 7'^. In the heavy baryon formulation this corresponds to replacing 2S!^ with v^ between spinor 
fields. 

This affects the determination of the constant Cpp ■ The partial width of a negative parity pentaquark is |23| 

r_(pa^5C^e)^(j^ace)2C_PSms^2^^^ (30) 

2t:j'^ mp 

The width for negative parity differs from the width for positive parity only in the replacement of fc^ by iJ^fc^r. For 
the total width of the 6, we find 

r_(e) == (966MeV)C2pB. (31) 

Thus, for negative parity, recent experiments suggest that C-ps is at most |, and the reanalysis of old KN scattering 
data suggests that C-pp is smaller than ■^. In other words, the small coupling expansion is even more reliable for a 
pentaquark of negative parity. Of course, we could look at this another way. It would be surprising for a pentaquark 
of very small width to have negative parity, as it requires an extremely small coupling and a corresponding dynamical 
mechanism. 

At NLO, the quark mass enters the pentaquark Lagrangian: 

^f = bn^PvM^Pv + '^P^(Tr M^P^P, + ... (32) 

Here Bq is the constant appearing in the meson Lagrangian in eq. H20() . The factor of -^ has been included so 

that NLO mass shifts are in terms of an order one coupling and -j-j^ rather than mq. An alternate definition of 
the coupling constants b-u and cp would be to leave out this factor. The dots in eq. (|32|l refer to terms that do not 
contribute to the self-energy at NLO. 

The quark mass contributes to the pentaquark self-energy at NLO as shown in figure^. We work in the isospin 
limit, with rn^ — rrid — m. Isospin breaking, by quark mass differences and by electromagnetism, is negligible 



numerically, until a higher order in chiral perturbation theory. In the isospin limit, the quark masses are related to 
the squared pion masses by 



2 - .4 2 



m^ = Bo{-ms + -to), to^ = Boirris + to), m^ = 2i?om. (33) 



^3 3 

For the pentaquark self-energy diagram in figure H^, we find 
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where the index a runs over pentaquarks. Summing powers of the self-energy in the usual way shifts the mass in the 



propagator, nia — rnp -f Sa . The constants a^ and a° in eq. H34|) are 
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-2hn - 2ap a^ ^ bn - <7p (35) 



C^K = -o^W ^ 2fTp 



a^ = --^^w - 2crp 
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There are four isospin multiplet masses and two parameters. Thus, there are two mass relations analogous to Gell- 
Mann's equal spacing rule for isospin multiplcts in the decuplet, 

TO-H — TOs == rn-^ — niN ~ ttin — rriQ. (36) 

Using the mass of the candidate S~~ seen in Q at 1860 MeV and the relation 

m^-me = bn { ^Y^f" ) , (37) 



we find bu = 1.17. Models for the pentaquark anti-decuplet [la, lla suggest that the constant b-j-c of order 1. The 
ap term in the Lagrangian only contributes to an overall shift to the masses at this order; however, it cannot be 
eliminated entirely since M^ also contributes vertices with pions. 

The mass relations obtained here differ considerably from those in [l5j because we do not consider mixing between 
the anti-decuplet and a degenerate octet. In particular, without mixing, the O is the lightest member of the anti- 
decuplet at 1540 MeV. The Roper resonance at 1440 MeV cannot be accommodated in the anti-decuplet. We will 
take up the inclusion of a degenerate octet in chiral perturbation theory elsewhere. 

III. POWER COUNTING 

To calculate to higher orders consistently, we need to establish the power counting for this theory. The coupling 
CpB was estimated to be j or smaller. Since it is about the same size as t^j "^6 will expand jointly in the chiral 
order and the order in the coupling constant. The order of a generic diagram is 

D^AL-2h-lH + J2k{N^ + Nli)+Nc. (38) 

fc 

The last term in this equation simply accounts for each of the Nc instances of a small coupling. It includes the 
small pentaquark-pion-nucleon coupling CpB and possibly other couplings expected to be small by our dynamical 
assumption, such as be in eq. (|42|l . The remaining terms in the equation account for the chiral order of the diagram. 
With a mass-independent renormalization scheme, the only dimensional quantities appearing are the small momenta. 
Thus, it is sufficient to determine the factors of momentum appearing in the diagram. There is a factor of 4 for each 
loop from d'^p, a factor of —2 for each internal meson line from the propagator, and a factor of —1 for each internal 
heavy particle line from the propagator. There is a factor of k for each of the N^ vertices from the fc**^ order meson 
Lagrangian and for each of the N^ vertices from the fc*^ order heavy particle Lagrangian. These terms account for 
the chiral order of the diagram. 



In addition, every diagram satisfies the identities, 



2Ih + Eh=J2 2^^- 



(39) 
(40) 



Here Eh is the number of external heavy particle lines. The first of these equations is an analogue of Euler's formula 
for a Feynman diagram; the second equation simply counts the "ends" of heavy particle lines in two different ways. 
Combining these three equations yields the power counting rule, 



1 



D^2L + 2--EH+y2[ik- l)iV^ + (fc - 2)N^] + Nc- 



(41) 



If all the heavy particles are nucleons and there are no instances of the small coupling, this reduces to the power 
counting rule for nucleon chiral perturbation theory. Since each term in eq. (|41|) is nonnegative, the diagrams 
contributing at a given order can be enumerated systematically. 

For example, the quark mass insertion in figure ^^ contributes to the mass splitting of pentaquarks at second order 
as already discussed in cq. 134|) . At third order, loop diagrams arise; however, they can only contain vertices from the 
first order Lagrangian in eq. H21() . There are two such diagrams with third chiral order, as in figure [2| The nucleon 
on the internal line in the second diagram suppresses it by two factors of the small coupling CpB relative to the first. 
The calculation of the diagrams is taken up in the next section. 



IV. MASSES AT NNLO 



The power counting developed in the previous section determines what diagrams could contribute to the mass shift 
at NNLO. There are one-loop diagrams with vertices from C^^' , tree diagrams with vertices from C^^' and at most 
two vertices from C^'^' , and tree diagrams with vertices from C^^' and at most one vertex from C^^' . Fortunately, 
the actual number of diagrams that contribute is small. Vertices in C^"^^ and £'^' that involve derivatives acting 
on pentaquarks fields contribute to the kinetic energy but do not affect the mass shift, i.e. the self-energy at zero 
momentum. In other words, the relevant terms in the higher order Lagrangians involve only the quark mass. Since 
Mq is second order, terms with factors of Mq and no derivatives appear only in the even order Lagrangians. The 



relevant terms in Cp appear in eq. H32|l 



(2) 

and the relevant terms in Cpg are 
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(42) 



It would seem then that there are two diagrams contributing to the mass shift at NNLO. There is the diagram in 
figure and the one-loop diagram in figure l^t- (The analogue of figure |3| with the internal baryon line replaced by 
a pentaquark is not one-particle-irreducible; its contribution is already included in the NLO mass shift.) The first of 
these is proportional to 6^. Note the similarity of the be term in eq. (|42() with the CpB term in eq. (|22(l . In fact, quark 
models suggest that the constant be should be small like CpB- Since the diagram in figure 13 is further suppressed by 
this small constant, it is higher order. That leaves one one-loop diagram contributing to the mass shift at NNLO. 

The diagram in figure |2Ji evaluates to 



-4^^^) = ^^ 



P.S'^StP.Y.GafKl^{-v-k) 
f 



(43) 



FIG. 2: Self-energy corrections appear at third order in the chiral expansion, 
nucleons, and dashed lines are pions. 



Double lines are pentaquarks, single lines are 



b) 



FIG. 3: Two insertions of the quark mass term in L 



(2) 



where a and h label pentaquarks and / is summed over mesons. The object P^ in this equation is just a spinor, while 
Gaf has all the flavor dependence, 

Gaf = (P" ^^■'^Tr^fe'ps y/)(P^'""^7r/'™P\p,), (44) 

where all repeated indices other than a and / are summed over. The function K comes from the loop integral, 

Calculating the loop integrals at zero external momentum, we find 

-^5]i3)(0) = '^T^S^P^^Y.Gr (£-?) . (46) 



/ 



Thus, up to NNLO the mass shift is 



m. = mp + <(£^)m^ + <(£y)m. - /3,"(£j)^m, - H'kC^f^K - ^l(^?^^ (47) 

where cA are the NLO shifts in eq. Ij^SI) and /?? — 27r7ipG/°°. The similarity between the forms of the NLO and 
NNLO terms is a bit misleading. The NLO terms are often expressed in terms of the quark mass. We have adopted 
a convention for the quark mass term in the pentaquark Lagrangian that allows us to use the LO relations for the 
squared pion masses in eq. H33|) to write the NLO terms in eq. 147|l in terms of the pion masses. The pion masses 
appearing in the NNLO terms on the other hand are physical, coming from pion loops. The same is true for the 
factors of 47r/; those in the NLO terms are conventional while those in the NNLO terms are physical, coming from 
pion loops. Of course, it's possible to turn this around and write both NLO and NNLO terms in terms the quark 
masses. As an example, we do this for the O mass, 



rrie = "^p — 2(6-^ + (jp)nis + 2(6-^ — 2ap)m 






4 2 3 3 

{—nig H — ?fi)2 + (jTis + m)2 



(48) 



where the tildes reflect an alternative definition of the coupling constants. The calculation of flavor factors such as in 
eq. (|44|1 for loop diagrams is greatly simplified using Mathematica. The result for /3? is 

(3^-Y'^I P%-y^% P?-0 (49) 

oN ""0^2 r)N n^n-i"^ rN ■"'0^2 

Pr, - l^Hp Pk = ^T^rip Pn ^ l^rLp 



At = ■k'^p Pk = -^'^p Pt = — ^i 



In the limit of massless up and down quarks, the numbers in eq. (|49|) reduce to an inverted form of the corresponding 
numbers for the decuplet |20j up to a constant factor coming from the spins. 

To summarize, at leading order the pentaquark masses are equal. The average mass of the anti-decuplet is mp and 
the spacing between isospin multiplets is 0. At NLO, the h-j-i quark mass term in the Lagrangian produces an equal 
spacing between isospin multiplets. The mass difference rriN — ttiq and its two analogues are nonzero and equal to 
5^^o where 

SNLO^^^^^!A_^^j. (50) 



A measurement of the spacing and eq. H50(l determines the constant b-n at NLO. As mentioned earher, using the mass 
of the candidate S (1860) seen in [g, we find b-j-^ = —1.17 at NLO. In addition, both the by^ term and the ap term 
shift the average mass of the anti-decuplet to 

-fi° = mp- {bn + 3^P) (^^%^) • (51) 

The data on masses are insufficient to determine mp and ap independently. Using eq. (|35|1 . we find 

mp-ap{ ^ j == 1860 MeV (52) 

At NNLO, the pion- induced self-energy shifts the masses. There is a shift in the average mass, a contribution to the 
average spacing between isospin multiplets, and a correction to equal spacing. The average mass of the anti-decuplet 
at NNLO is 

., 'ml + 4:171% -|-3m?.\ , ^ 

Although the NNLO shift is suppressed by one chiral order or about i it turns out that the self-energy integral 
contributes an extra 2tt. Thus, the NNLO shift in the average mass could be comparable to the NLO shift. This 
simply means that the NNLO fit for the combination of mp and ap in eq. (|52|l could differ from the NLO fit more 
than expected. 

The average spacing between isospin multiplets at NNLO is 

gNNLO ^ SNLO + ^2^ /'M±!<_^'\ (54) 

The NNLO contribution to the average spacing could also be comparable to the NLO contribution. Thus, the NNLO 
fit for bn could differ from the NLO fit more than expected. 

At NNLO, the difference (ms — itin) — {rriN — ttiq) and its two analogues are nonzero and equal to g^^^o -^^j^gj-g 

MNLO__ ^2 ( iml-4m\ + ml \ 

= 7i:p (1.31 MeV). (56) 

Unlike the NNLO effects in eq. (j53|l and eq. (|54|l . c^^^o jg suppressed significantly relative to NLO effects like S^^'~' . 
In principle, a measurement of e^^^o dgtermincs Hp, which can then be used to determine b-j-c at NNLO via eq. H54|l. 
Current data are insufficient to carry through this analysis. As mentioned earlier, based on quark models we expect 
that Hp is of order 1. In this case, the correction to equal spacing is roughly IMeV. In other words, we expect 
the corrections to the equal spacing rule to be negligible and difhcult to measure experimentally. It is possible that 
such a cancellation does not occur at NNNLO. Dimensional analysis and a renormalization scale of /x = Airf suggest 
corrections of about 15 MeV. 

Thus, the problem of determining Tip, and of determining b^. at NNLO, remains. In addition, unlike CpB, "Hp 
cannot be determined using decays because decays within the anti-decuplet are kinematically forbidden. One possi- 
bility that needs to be explored further is the determination of Tip and bn -\- 'i.ap on the lattice using the behavior of 
the pentaquark mass as a function of the light quark mass in eq. I|48|l . Lattice work on pentaquarks is well underway. 
For example, Csikor et. al. identify a negative parity candidate for 0+ [lj| . 

The large but benign NNLO effects described above were originally studied in the nucleon and decuplet sector in 
|2Clj| . In that paper, it was also pointed out that, in the limit of massless up and down quarks, the difference g^^^c* 
contains a suppression factor of (-7= — 1). For example, 

NNLO _ Tj2 /^ 2 \ m^ 



which is about 20 MeV for Tip of order 1. In 21], flavor 27 contributions to the mass, such as those coming from pion 
loops, were studied in greater generality. It was noted that the suppression of flavor 27 contributions is even greater 
with nonzero up and down quark masses, as seen in eq. H55|l above. 
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At one higher order, there are contributions from one-loop diagrams with quark mass insertions on the internal line, 
from wavefunction renormalization of the external lines in the tree level diagram, and from terms in the fourth order 

4 2 

Lagrangian with two factors of the quark mass. The first two of these contribute logarithmic terms like ,^™-?^j In (^t^) 
to the mass shift. One- loop diagrams and quark mass insertion diagrams with two factors of a small coupling such as 
in figure |2l3 and figure |31 appear at one higher order yet. 

V. NEGATIVE PARITY 

The calculation of the mass shift can be easily repeated for pentaquarks of negative parity. As discussed earlier, the 
Lagrangian Cp is unchanged under a change in parity of the pentaquark field because the field appears quadratically 
in each term. The Lagrangian CpB is altered, but this affects only the determination of the constant Cps- 

The mass shift calculation proceeds exactly as in the positive parity case. Only the self-energy diagrams with 
pentaquarks coupling to non-exotic baryons, in figures l^b and 13 are affected. As pointed out in the previous section, 
these diagrams are fifth order. Thus, the mass shifts are not sensitive to parity even at NNNLO, which includes 
leading logarithms. 

This suggests that negative and positive parity states should behave similarly under chiral extrapolation. In 
particular, the ordering between these states should remain unchanged under chiral extrapolation. This result holds 
even if a pentaquark octet is included, as long as it is of the same parity as the anti-decuplet. 

VI. CONCLUSION 

Heavy baryon chiral perturbation provides a framework to calculate the mass splitting in the pentaquark anti- 
decuplet systematically up to NNLO. The calculation could be complicated by baryons in the nearby energy range. 
However, an expansion in the coupling of pentaquarks to non-exotic baryons simplifies calculations. With this coupling 
expansion it turns out that the mass splitting is independent of parity even at NNNLO. Corrections to the equal spacing 
rule at NNLO are negligible. Our results for the pentaquark mass splitting as a function of the pion masses should be 
helpful in chiral extrapolation of lattice results. In addition, lattice methods could help determine the two coupling 
constants in the chiral Lagrangian. In future work, it will be important to consider the effects of the other pentaquark 
multiplets predicted by spin-flavor symmetry in diquark models, uncorrelated quark models, or large N . 
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